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We present recursive multiport schemes for implementing quantum Fourier transforms and the 
inversion step in Grover’s algorithm on an integrated linear optics device. In particular, each 
scheme shows how to execute a quantum operation on 2d modes using a pair of circuits for the 
same operation on d modes. The circuits operate on path-encoded qudits and realize d-dimensional 
unitary transformations on these states using linear optical networks with O (d^) optical elements. 

To evaluate the schemes against realistic errors, we ran simulations of proof-of-principle experiments 
using a simple fabrication model of silicon-based photonic integrated devices that employ directional 
couplers and thermo-optic modulators for beam splitters and phase shifters, respectively. We find 
that high-fidelity performance is achievable with our multiport circuits for 2-qubit and 3-qubit 
quantum Fourier transforms, and for quantum search on four-item and eight-item databases. 

PACS numbers: 03.67.Ac, 42.50.Ex 


I. INTRODUCTION 

Linear optics with single photon sources and detec¬ 
tors provides a promising candidate for establishing effi¬ 
cient and scalable universal quantum computation mm- 
Some specific advantages of optical implementations are 
the robustness of photons against decoherence and ultra¬ 
fast optical processing. Meanwhile, current limitations 
include low efficiencies in photon creation and measure¬ 
ment, and the difficulty of storing photonic quantum in¬ 
formation in a quantum memory. However, the main 
technical challenge with the scalability of linear optical 
devices is the considerable overhead necessary for realiz¬ 
ing multi-qubit gates. Moreover, maintaining phase sta¬ 
bility between optical modes often requires active check¬ 
ing and calibration—a task that clearly becomes more 
demanding the larger the circuit gets. Both of these is¬ 
sues are addressed by integrated photonic technology. 

A photonic integrated circuit (PIC) is a multiport de¬ 
vice consisting of an integrated system of optical elements 
embedded onto a single chip using a waveguide architec¬ 
ture US]. Because PICs are compact and designed to 
have inherent phase stability, they offer the potential for a 
truly scalable optical quantum computer. They are also 
fully compatible with electronic devices and fiber optic 
systems, which can lead to increased functionality. 

There has been much progress made in PIC-based ap¬ 
proaches to quantum optics: recent experiments have 
demonstrated linear optical quantum gates 016!, multi¬ 
photon entanglement 00 boson sampling experiments 
0{n], quantum walks in optical arrays [SHni, and 
simulation of quantum systems |151 I16j . More gener¬ 
ally, PICs offer a natural platform for conducting ex¬ 
periments on quantum systems with higher-dimensional 
Hilbert spaces [niiiEi- 
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In this paper, we describe recursive multiport schemes 
for carrying out quantum operations on a PIC. By re¬ 
cursive we mean that two copies of the d-dimensional 
circuit are used to construct the 2d-dimensional version. 
Formally, a multiport circuit represents a decomposition 
of a unitary transformation into a network of single¬ 
qubit gates acting on adjacent modes, which on a PIC 
is mapped onto a sequence of beam splitters and phase 
shifters. It is crucial to note that in our schemes, the 
photons represent multi-rail qudits, not dual-rail qubits. 
This is an important distinction since existing methods 
for performing entangling gates on dual-rail qubits are 
not very scalable. However, our circuits act on single 
qudits, therefore, such entangling gates are not needed. 

We must emphasize, though, that linear optical imple¬ 
mentations on single photonic qudits are inherently un¬ 
scalable for quantum computation since they require an 
exponential number of optical modes. The only known 
way for performing scalable quantum computation with 
linear optics requires multi-photon quantum interference 
and and measurement-based nonlinearity, such as in the 
KLM proposal 0 . Nevertheless, there are several con¬ 
texts in which our recursive circuits may find suitable 
application, for instance, within a larger architecture of 
entangled qudits HU or when using QFT as a verifica¬ 
tion tool in boson sampling [5D|, the latter having been 
experimentally demonstrated by Carolan, et al. |2Ij . 

Here, we consider circuits for two important families 
of unitary transformations: quantum Fourier transform, 
which is an important subroutine in many quantum algo¬ 
rithms, is covered in Section jnj and inversion about the 
mean, which is a key step for the iterations in Grover’s 
algorithm, is covered in Section |HI[ In both cases, we 
provide a recipe for constructing a linear optical net¬ 
work that uses a total of O (^cP) beam splitters and phase 
shifters for realizing a quantum operation on d modes. 
The basis for these constructions are efficient matrix fac¬ 
torizations of the 2c?-dimensional unitary operators into 
a product of a matrix consisting of two copies of a d- 
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dimensional version of the same unitary operator and 
some sparse matrices with 2 x 2 blocks on the diago¬ 
nal. Thus, our schemes provide a systematic way of con¬ 
structing larger circuits from smaller ones, which is quite 
beneficial for scaling these quantum operations to higher 
dimensions. 

To demonstrate the practical viability of our multiport 
scheme, we also conducted simulations of experiments on 
the circuits for quantum Fourier transform and Grover’s 
algorithm using a fabrication model that incorporates 
realistic errors in the beam splitters and phase shifters 
based on wafer-scale testing data on PIC components 


The simulation results are discussed in Section IV 


II. QUANTUM FOURIER TRANSFORMS 



FIG. 1: (Color online) The optical multiport circuit for 
2-qubit QFT F 4 . 

For simplicity, B(i,j) denotes an equal beam split¬ 
ter (e = 1 / 2 ) between modes i and j. 

2. S{i,j) refers to a swap operation between modes i 
and j. With the convention above, this is the same 
as having a beam splitter with e = 0 . 


Many known quantum algorithms that exhibit expo¬ 
nential speedup over their classical counterparts make use 
of quantum Fourier transform (QFT), which describes a 
discrete Fourier transform on quantum mechanical am¬ 
plitudes [23) . Typically, the speedups come from per¬ 
forming quantum phase estimation, which involves find¬ 
ing approximate eigenvalues of a unitary operator, and it 
involves an inverse QFT. It is especially helpful in solving 
interesting problems like prime factorization and discrete 
logarithms 

There have been some recent experiments that real¬ 
ize QFT with optical multiport circuits, in particular, 
the four-mode version by Laing et al. |25j . which was 
performed with path-and-polarization encoded states in 
bulk optics, and a six-mode version used in the study 
of inequivalent classes of complex Hadamard matrices by 
Garolan, et al. m- These examples provide evidence 
for the considerable interest in developing practical lin¬ 
ear optical methods for implementing QFT. 

In this section, we describe a recursive multiport 
scheme for realizing QFT with integrated linear optics. 
The key ingredient for the matrix decomposition involved 
is the unitary operator given by the direct sum of two of 
Fourier matrices. 

In order to facilitate the description of multiport cir¬ 
cuits, we shall use the following special notation through¬ 
out this paper. A multiport circuit C with d optical 
modes is called (7(1, 2,..., d). The modes are labeled 1 
to d from top to bottom. Optical elements are organized 
from left to right according to the sequence they appear 
in the circuit. Those that can be performed in parallel 
are enclosed in square brackets. 

We will denote the optical elements as follows: 

1 . refers to a beam splitter with reflectivity e 

acting on modes i and j. Our convention here is 
to choose the overall phases so that in terms of the 
mode operators we have 


( 4out\ = ^ 

Vslout/ 


( 1 ) 


3. Pe{i) refers to a phase shifter on mode i, which 
means the amplitude for mode i is multiplied by 


Using our multiport circuit notation, the single-qubit 
QFT, or the Hadamard gate, is just ^ 2 ( 1 , 2) = i?(l,2). 
For 2-qubit QFT, the multiport circuit is given by 

F4(1, 2,3,4) = [5(2,3)] [B{1, 2) 5(3,4)] [P^ (4)] 

[5(2,3)] [5(1,2) 5(3,4)] [5(2,3)]. (2) 

The circuit for F 4 is illustrated in Fig. [l] 

Let d = 2’" be the number of modes for the m-qubit 
QFT multiport circuit. Using F 4 we can construct the 
3-qubit QFT Fg as follows: 


58(1,2,..., 8) = [5(2,3) 5(4,5) 5(6,7)] [5(3,4) 5(5,6)] 
[5(4,5)] [ 54 ( 1 , 2 , 3 , 4 ) 54 ( 5 , 6 , 7 , 8 )] 


55(6)55(7)53,(8) 


]5(4,5)] [5(3,4) 5(5,6)] 


[5(2,3) 5(4,5) 5(6,7)] 

[ 5 ( 1 , 2 ) 5 ( 3 , 4 ) 5 ( 5 , 6 ) 5 ( 7 , 8 )] 

[5(2,3) 5(4,5) 5(6,7)] [5(3,4) 5(5,6)] 

[5(4,5)] (3) 


Fig. [2] shows a diagram for how circuits for F 4 is used 
for performing the 3-qubit QFT Fg. 

To describe how to build F 2 d from 5^, it is convenient 
to define the shuffle operation E that performs the fol¬ 
lowing permutation on 2 d modes: 


( 1 , 2 ,..., 2 d)^ 

(1, d 4- 1,2, d -j- 2,..., A:, d 4- A:,..., d, 2d). (4) 


Let E ^ denote the inverse permutation. 

Observe that E(l, 2,..., 2d) can be realized with a PIG 

d(d-l) , 

using -^-- 2 —- swap gates: 

E(l, 2,..., 2d) = [5(d, d -f l)][5(d - 1, d) 5(d +l,d + 2)] 
[5(d - 2, d - 1) 5(d, d-4l)5(d-4 2,d-4 3)]--- 
]5(2,3)5(4,5)---5(71-2,11-1)]. 


( 5 ) 
















3 



FIG. 2: (Color online) Implementing 3-qubit QFT using a pair of multiport circuits for 2-qubit QFT. 


Thus, 

F2d(l,2,...,2d)= [E-i(l,2,...,2d)] 

[Fd{l,2,...,d)Fd{d+l,d + 2,...,2d)] 

P^id + 2)---Pk^id + k + l)--- P (^-i). {2d) 
[S(l,2,...,2d)] 

[B{1,2) B{Z,F) ■ ■ ■ B{2d - \,2d)] 
[E-i(l,2,...,2d)] (6) 

Let N{Fd) denote the number of optical elements 
needed for the d-dimensional QFT multiport circuit. 
With N{F 2 ) = 1 and N{F 4 ) = 8, we get the following 
recursive formula for the circuit size: 

lV(i^2d) = 2iV(Fd) + ^^^^^+d+(d-l). (7) 


The first term in the sum on the right-hand-side corre¬ 
sponds to a pair of QFT circuits on d modes, the second 
term counts the total number of swap gates, the third 
term counts the equal beam splitters, and the last term 
counts the phase shifters. Expressed as a function of d, 
we have 


N{Fd) 


3d2 -b d (log2 d - 7) 

4 ^ ^ 


( 8 ) 


so the total number of elements in the multiport circuit 
is quadratic in the number of modes, as expected. 

The recursive scheme works not just for qubits but for 
any QFT circuit with an even number of modes. More 
precisely, given an initial QFT circuit F^ the scheme 
can be used to implement QFT for d = n2^, where k is 
a positive integer. 

The scheme is essentially equivalent to a PIC transla¬ 
tion of an important matrix factorization of Fourier ma¬ 
trices, initially discovered by Gauss [55] and a precursor 
to algorithms that perform fast Fourier transform m- 


F2d — 


1 

71 






( 9 ) 


where F^ is the d-dimensional Fourier matrix, 1^ is the d- 
dimensional identity matrix, Dd = diag (l,a;,... ,oj'^') is 


a diagonal matrix with to = and P is the 2d x 2d 

permutation matrix that maps the column vector v = 
(vi,V 2 , ■ ■ ■,V 2 d)'^ into 

Pv = (ui,U 3 , . . . ,V2d-l,V2,V4:, ■ ■ . , U 2 d)^ , (10) 

that is, it shuffles components of v such that the first 
half involves components with odd indices and the second 
part involves the even ones. Note that P is actually the 
matrix representing the permutation 

It is important to note that linear optical implementa¬ 
tions of QFT have been previously explored by Torma, 
et al. [55], and Barak and Ben-Aryeh |53|. Torma, et al. 
examine the sufficient number of beam splitters needed 
for totally symmetric mode couplers, of which the dis¬ 
crete Fourier transform is a special case. Their calcula¬ 
tions show that for a d-mode circuit, (dlog 2 d )/2 beam 
splitters are sufficient. On the other hand, Barak and 
Ben-Aryeh describe a particular linear optical scheme for 
QFT based on the Gooley-Tukey algorithm. 

A crucial difference between our scheme and these 
other approaches is that we are restricted to beam split¬ 
ters that operate only on adjacent modes, since this is 
a limitation on PIGs. This is why our scheme generally 
requires more beam splitters, in order to perform those 
additional swap operations. Both of the previous schemes 
were designed with bulk optics in mind, so they do not 
consider this restriction. 

For QFT, all three approaches are formally equivalent 
through Eq. [^ In fact, the same number of equal beam 
splitters is used in all schemes; it is the number of phase 
shifters that vary. 

Torma, et al. provide a general formula for the matrix 
factorization they used but it does not completely specify 
where phase shifts are actually needed, since they are 
mostly concerned with counting the beam splitters. 

Barak and Ben-Aryeh describe a specific Gooley-Tukey 
factorization of d-dimensional Fourier matrices into log 2 d 
unitary operations that contain d/2 pairs of beam split¬ 
ters and phase shifters, i.e., the phase shifts are always 
coupled to a beam splitter. We note that a recent exper¬ 
iment by Crespi, et al. |30| for testing the quantum sup¬ 
pression law |20] implements the Barak and Ben-Aryeh 
circuits for four- and eight-mode QFT on a 3-D PIC. 
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In contrast, our scheme employs d—3 phase shifters for 
d-mode QFT, so we achieve a modest savings on phase 
shifters. For example, the 3-qubit QFT circuit of Barak 
and Ben-Aryeh uses 12 phase shifters but ours uses only 
5. It may be worth mentioning that if one actually im¬ 
plements our QFT circuit with bulk optics, or a 3-D PIC 
such as in Ref. [SD], where the swap gates do not in¬ 
volve beam splitters, our circuit is slightly more efficient 
because it requires fewer phase shifters. 


III. THE INVERSION STEP IN GROVER’S 
ALGORITHM 

Grover’s algorithm |31j describes a quantum algo¬ 
rithm for searching an unsorted database of d items us¬ 
ing 0{Vd) calls to an oracle, which offers a quadratic 
speedup over known classical methods. In the most ba¬ 
sic scenario, we have a database with d items and we are 
supplied with a quantum oracle that can mark the solu¬ 
tion to the search problem by shifting the phase of the 
solution’s register. The goal of the algorithm is to find 
the solution using the smallest number of queries. 

To start, we prepare the equal superposition state 

( 11 ) 

x—1 

where d = 2 "* is the number of basis states |a:) for m 
qubits. 

Grover’s algorithm is then characterized by the re¬ 
peated use of a quantum subroutine known as the Grover 
operator 


G=(2|^)(V^|-/)0, 

( 12 ) 

where O is a query to the oracle and 


W = 2\iP){iIj\- I 

(13) 


is often called the inversion about the mean. We shall re¬ 
fer to W as Grover inversion. In this section, we describe 
a recursive multiport scheme for implementing Grover in¬ 
version on a PIG. In particular, the multiport circuit for 
searching d items is utilized as a building block to the 
circuit for searching 2 d items. 

If we consider a search problem with a unique solution, 
the oracle can be realized by a single 7r-phase shift on the 
optical mode corresponding to the item to be marked. 

To describe the main result, it is helpful to first con¬ 
sider the unitary transformation Vd, which involves a 
relatively simple network of equal beam splitters on d 
modes. Using the multiport notation presented in Sec¬ 
tion]^ we have 

V4(1,2,3,4) = [R(1,2)B(3,4)] [^(2,3)] 

[R(1,2) 5(3,4)] [5(2,3)]. (14) 

The multiport circuit for V 4 is depicted in Fig. 



FIG. 3: (Golor online) The optical multiport circuit for 
the unitary V 4 . 


From V 4 , we can build any other circuit Vd for d = 2^ 
modes according to the rule 

V2d{l, 2,..., 2d) = [Vdih ..., d) 14(d + 1,..., 2d)] 

[E(l, 2,..., 2d)][5(l, 2) 5(3,4) • • • 5(2d - 1, 2d)] 
[E-i(l,2,...,2d)] , (15) 

where S is the same shuffle operator defined in Eq. (§. 
As an example, the circuit for Vs is shown in Fig. . 

Let Wd denote Grover inversion on d modes. First let 
us consider W 4 , which is given by 


IU4(1, 2,3,4) = [5(1,2) 5(3,4)] [5(1,2) 5(3,4)] [5(2,3)] 
[5(1,2)] [5(2,3)] [5(1,2) 5(3,4)]. (16) 

Fig.[^shows how W 4 can be implemented with integrated 
optics. 

In the general case, it is convenient to define the uni¬ 
tary operation $ on 2 d modes that exchanges the photon 
amplitude in mode 1 and mode d -I- 1 through the follow¬ 
ing network of swap gates on neighboring modes: 


$(1,2,..., 2d) = [5(d,d-bl)] 

[5(1,2) 5(3,4) •••5(d-l,d)] 

[5(2,3) 5(4,5)---5(d-2,d- I)]--- 




^' 1-^4 


s\i 




[5(1,2) 5(3,4) • • • 5(d - 1, d)][5(d, d + 1)]. (17) 


It is worth mentioning that $(1,2..., 2d) employs a total 

,2 j 

of ^ -I -1 -I -1 swap gates. 

The general multiport circuit for Grover inversion W 2 d 
employs both Wd and Vd^ and is given by 


W2d{l, 2,..., 2d) = [Wd{l, ... d) Wd{d + 1,..., 2d)] 
[Vdil, ... d) Ud(d + 1,..., 2d)] [$(1, 2,..., 2d)] 
[Urf(l,...d)14(d+l,...,2d)]. (18) 

Fig. |6] illustrates how the rule is used for building the 
circuit for Ws- 

Gounting the number of optical elements used in W 2 d, 
first for the unitary Vd, we have 


N{V2d) = 2N{Vd)+d{d-l) + d, 


(19) 






















5 



FIG. 4: (Color online) Implementing the unitary Vg using two multiport circuits for V4. 



FIG. 5: (Color online) The optical multiport circuit for 
Grover inversion W 4 . 


and Q is the 2d x 2d permutation matrix that exchanges 
the first and {d + l)th entry of a column vector. Unlike 
in the QFT circuit, the construction for Grover inversion 
is known only to work when d = 2^, because we do not 
know of any natural counterpart to Vd when d is not a 
power of two. 


IV. MULTIPORT CIRCUIT SIMULATIONS 


where the first term refers to a pair of Vd circuits, the 
second term refers to the elements in S and and the 
last term refers to a set of parallel equal beam splitters on 
neighboring modes. Solving the formula with N (V 4 ) = 6, 
we obtain 


N{Vd) = ^^^^. ( 20 ) 

Now for Grover inversion, we have 
N {W 2 d) = 2N (Wd) + 47V {Vd) +d^+d+l, (21) 


where N{W 4 ) = 9. Plugging in N (Vd) and solving the 
relation yields 


iv(w,) = 57!:i4<5hlU±i)_i. 

8 


( 22 ) 


The recursive scheme directly implies that the Wd has 
a matrix decomposition given by 


W2d 


(Vd 0 \ / Ud 0 \ / lUd 0 \ 

1^0 Vd)^\Q Wd)' 




(23) 

(24) 


where Id is the d-dimensional identity matrix. 


IU2 = 



1 A 





(25) 


In realistic linear optical systems, optical elements ex¬ 
perience photon losses, optical modes suffer from relative 
phase mismatches, and fabrication defects lead to errors 
in the splitting ratios of beam splitters. To account for 
such device imperfections, we follow the example of Ref. 
[22] and consider a simple model for silicon-based PICs 
that use directional couplers [35] for beam splitters and 
thermo-optic phase modulators [33] for phase shifters. 

In this section, we simulate experiments on our mul¬ 
tiport circuits and assess their performance under a fab¬ 
rication model that focuses on two primary sources of 
errors: (i) incorrect reflectivities in beam splitters and 
(ii) absorption losses in the phase shifters. 

For directional couplers, the error in the splitting ratio 
is due to imperfections in the dimensions of the coupled 
waveguides. In our model, we assign reflectivities to the 
beam splitters in the multiport circuits according to a 
Gaussian distribution, with mean 0.5 and standard devi¬ 
ation 0.04, which agrees with the testing data on recently 
developed devices [33] . 

For thermo-optic modulators, the important source of 
error is the tree-carrier absorption in doped silicon ma¬ 
terial, leading to propagation loss. Typically, the ab¬ 
sorption process is modeled as a unitary operation by 
introducing a beam splitter between the lossy mode and 
an ancillary one, and whose transmissivity represents the 
photon loss rate for the device. Fortunately, the matrix 
representing a phase shifter is always block diagonal in 
the circuit. It is therefore sufficient to apply a scaling fac¬ 
tor y/\ — 7 on the lossy mode, where 7 is the absorptivity. 
In our model, we assign 7 values to the phase shifters ac¬ 
cording to a rectified Gaussian distribution, with mean 
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FIG. 6: (Color online) Implementing Grover inversion Wg using pairs of multiport circuits for W 4 and V4. 


0.05 and standard deviation 0.025, consistent with the 
reported loss rates in the latest studies |55] . 

Our scheme also utilizes a significant number of swap 
gates, each of which can be realized by a beam splitter 
with vanishing reflectivity. However, it may be possi¬ 
ble to implement the swap more efficiently since there is 
no required interaction between the modes. As such, we 
chose to model them as slightly better performing beam 
splitters, whose reflectivities are drawn from a rectified 
Gaussian distribution of mean 0.02 and standard devia¬ 
tion 0.02. 

The first experiment we considered involves imple¬ 
menting the QFT circuit with randomly generated inputs 
of the form 


|(() 4 ) = (Z1,Z2,Z3,Z4)'^ /(4 

in the 2-qubit case and 

l^s) = (z:i,z2,z3,z4,z5,ze,zr,zs)^ /(a 

in the 3-qubit case where Zi = \nxi exp (27rj/i) and 
Q = X]i=i The parameters Xi and yi were all picked 
independently and uniformly at random from the interval 
( 0 , 1 ); 

It is known that a Haar-random d-dimensional pure 
state may be constructed by normalizing a vector of d 
independent, identically distributed complex Gaussian 
random variables |34j . Thus, the above procedure for 
generating input states is equivalent to sampling random 
pure quantum states from the uniform Haar measure, 
since Zi represents the Box-Muller transform [35) . which 
generates a pair of Gaussian random numbers from a pair 
of uniformly distributed ones. 

To evaluate the performance of our circuits, we calcu¬ 
lated the (squared) fidelity F(4',<i)) = | ('I'|$) p between 
the simulated output state |<i)) and the ideal output Idt) 
for each particular input. Because of the Haar-random 
sampling of input states, the mean fidelity we get pro¬ 
vides a good estimate of the average gate fidelity for the 
QFT circuit, a fairly common figure of merit for quantum 
gate experiments. 


After running 10^ trials, we obtained fidelities with a 
mean value of 0.94 and a standard deviation of 0.032 in 
the 2-qubit case, and a mean value of 0.86 and a standard 
deviation of 0.056 in the 3-qubit case. 

The second experiment involves running Grover’s algo¬ 
rithm on four-item and eight-item databases, each with 
a unique solution. In contrast with the first experiment, 
Grover’s algorithm specifically uses the equal superposi¬ 
tion state as its input. Thus, we have included the gates 
required for preparing this state in our simulation. 

The multiport circuit for a four-item Grover search is 
illustrated in Fig. The schematic for the eight-item 
version is displayed in Fig. where the dashed boxes 
refer to oracle queries realized by a 7r-phase shift on the 
appropriate mode, and the circuit Pg for preparing the 
input is shown separately in Fig. 

Observe that the oracle query and Grover inversion are 
repeated twice in the 8-item quantum search, which fol¬ 
lows from the fact that Grover’s algorithm prescribes do¬ 


ing 



iterations, which is 2 when d = 8 . In this case. 


Grover’s algorithm produces an output state that gener¬ 
ates the solution with high probability ~ 0.945), not 


deterministically. 


We performed 10^ trials for both searches, each with 
a randomly selected unique solution. For each trial we 
computed the fidelity of the simulated output with the 
ideal one. We obtained fidelities with a mean of 0.90 
and a standard deviation of 0.051 for a four-item Grover 
search, and a mean of 0.76 and a standard deviation of 
0.099 for an eight-item Grover search 

We must emphasize that what these results indicate 
is the feasibility of a classical simulation of Grover’s al¬ 
gorithm with our circuits. Nevertheless, like QFT, it is 
possible that the unitary Wd may find suitable applica¬ 
tion in some other truly quantum setting. 

There have been no linear optical experiments to date 
to which we can directly compare our simulation values; 
however, recent demonstrations on PIG devices have ob¬ 
tained fidelities of more than 0.94 mm- Gonsidering 
that our simulations were done with a rudimentary yet 
conservative error model, it is encouraging to see a com- 
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FIG. 7: (Color online) The optical multiport circuit for 
Grover search on four items. In the diagram, the second 
mode represents the solution to the search problem. 
The output of the algorithm is obtained by a final 
measurement in the standard basis. 
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FIG. 8: (Color online) The optical multiport circuit for 
Grover search on eight items. The box Pg corresponds 
to the circuit in Fig. In this diagram, the third mode 
represents the solution to the search problem. The 
output is obtained by a final measurement in the 
standard basis. 


parable level of performance with our multiport circuits. 


V. DISCUSSION 

Any d-dimensional unitary operation can be imple¬ 
mented with linear optics, by a triangular array of beam 
splitters on d modes, each accompanied by a phase 
shifter, and another d — 1 phase shifters [35]. Such a 
matrix decomposition requires d^ — 1 optical elements, 
exactly the number of real parameters needed to spec¬ 
ify a unitary in SU(d). While this specific arrangement 
might be preferred when trying to develop a fully recon- 
figurable linear optical processor, it is generally not the 
optimal choice for specific unitary families. If we wish to 
develop a quantum processor dedicated to a particular 
task, it may be possible, and quite desirable, to design a 
multiport circuit that uses fewer elements, since fabrica¬ 
tion defects in devices are known to severely hamper the 
quantum performance of PIGs. Furthermore, depending 
on how technological capabilities improve, photonic cir¬ 
cuits that consist just of equal beam splitters and a small 



FIG. 9: (Golor online) The optical multiport circuit Pg 
for preparing an equal superposition state for eight 
modes given an input photon in the first mode. 


Experiment 

Mean fidelity 

Std. deviation 

2-qubit QFT [8] 

0.944 

0.0319 

4-item Grover [14] 

0.904 

0.0507 

3-qubit QFT [41] 

0.861 

0.0559 

8-item Grover [112] 

0.762 

0.0990 


TABLE I: Simulation fidelities for 10^ trials. The total 
number of beam splitters and phase shifters in each 
circuit is indicated in the brackets. 


number of simple phase shifts might be more appropriate 
in certain contexts. 

The recursive scheme described in this paper illustrates 
that for unitary operations that serve as major compo¬ 
nents of Shor’s and Grover’s algorithm, we can assemble 
a multiport circuit from smaller versions of the same op¬ 
eration. The matrix decompositions for QFT and Grover 
inversion both have a relatively simple structure, which 
is made seemingly complicated only because of permuta¬ 
tions have to be implemented by a sequence of nearest- 
neighbor swap gates. 

We have previously discussed the optical circuit size 
for QFT and Grover inversion. For completeness, we can 
also remark on the optical circuit depth in comparison 
to the unitary matrix factorization of Reck, et al. Let 
D(U) denote the circuit depth for implementing U. In 
their case, the configuration of beam splitters and phase 
shifter is fixed for any arbitrary d-dimensional unitary so 
the circuit depth is 2{d — 1) . For the QFT circuits, we 
have 


D{Fd)=3{d-l)-2\og^d (26) 

while for the Grover inversion circuits, we have 

D{Wd) = 5d-2\og2d-(log^df-6. (27) 

Thus, there is a linear dependence on d in all cases, 
although our circuit depths are around a factor of 2 
worse. This is precisely due to the swap operations in 
our scheme, which is quite unavoidable because elemen¬ 
tary gates on a PIC are restricted to adjacent optical 
modes. 
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FIG. 10: (Color online) Histograms of simulation results for (a) 2-qubit QFT (0.949), (b) four-item Grover search 
(0.912), (c) 3-qubit QFT (0.870), and (d) eight-item Grover search (0.774), with bin sizes obtained using the 
Freedman-Diaconis rule . The median values are indicated in parentheses. 


The simulation results listed in Table U show that these 
circuits perform well even when realistic errors in the op¬ 
tical elements are taken into account. This is accom¬ 
plished even before we apply post-fabrication optimiza¬ 
tion techniques |22j that would significantly improve the 
overall performance. 

For a quick summary of the data, we include his¬ 
tograms of the fidelities in Fig. 10 where the bin size 
is determined by the Freedman-Diaconis rule. This gives 
a rough picture of how the fidelity values are distributed, 
in particular, how much are the values skewed to the left. 

Nonetheless, it seems clear that further advances in 
fabrication methods will be needed to scale to much 
larger circuits. From the results of both experiments, 
we can roughly estimate a 9% to 15% reduction in aver¬ 
age fidelity when the number of modes is doubled. This 
seems to be the case even though the number of elements 
in QFT increases fivefold while that of Grover’s algorithm 
increases eightfold. Still, this rate of decline implies the 
fidelity will go below acceptable levels rather quickly. It 
is previously known that errors of as high as 1% per ele¬ 
mentary gate will still allow fault-tolerant quantum com¬ 
puting m but in terms of the overhead resources needed 
to protect against failures, we will likely require less than 
0.1 % error per gate for a practical device. 

There are other long-standing issues with integrated 
photonic platforms. One is the problem of integrating 
single photon sources and detectors with a silicon-based 
PIG. Some recent progress in this area include waveguide- 
integrated semiconductor quantum dot sources [35] and 
superconducting nanowire detectors |39j . 


Another problem is that one typically needs to build 
customized chips for each experiment. Furthermore, 
settings may need to be reconfigured between different 
runs of the same experiment; for example, in an adap¬ 
tive experiment, several device settings might be up¬ 
dated according to the outcome of intermediate mea¬ 
surements. These latter concerns will be addressed by 
a reprogrammable photonic quantum processor Emu]. 


VI. CONCLUDING REMARKS 

Photonic approaches to quantum information process¬ 
ing has seen extensive growth over the last few years and 
linear optics has played a major role in many of the re¬ 
cent advances. In particular, optical implementations on 
integrated systems have vastly improved our ability to 
perform quantum experiments with single photons. 

In this paper we described recursive multiport circuits 
for quantum Fourier transforms and Grover inversion. 
The circuits are designed to be implemented on linear 
optical networks where qudits are encoded in the path 
traversed by single photons. Each circuit requires 0{cP) 
optical devices for implementing a d-dimensional unitary 
operation, which shows that they scale in a reasonable 
manner. We also demonstrated that these circuits can 
achieve high-fidelity performance in a practical setting by 
conducting simulations on a silicon-based PIC model that 
incorporates the faulty operation of real optical devices. 

Our recursive schemes take advantage of the natural 
modularity of photonic integrated circuits, allowing us 
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to construct more sophisticated multiport circuits using 
smaller versions of the same operation as building blocks. 
This may prove beneficial for practical implementations 
because it provides a straightforward approach to scaling 
operations into higher dimensions. In contrast, if one for 
example uses the triangular array of Reck, et al. [36] , the 
phases and reflectivities in the circuit will generally re¬ 
quire adjusting whenever more optical modes are added. 

Our results show the versatility of PICs as a platform 
for developing high-quality, scalable quantum informa¬ 
tion processors. It is our hope that this work will help 
promote new avenues of research into integrated photonic 


devices and stimulate technological improvements that 
will allow more practical applications. 
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